Method for estimating the transfer function of a multicarrier signal transmission channel and corresponding receiver

ABSTRACT

The invention concerns a method for estimating the transfer function of a multicarrier signal transmission channel consisting of a temporal succession of symbols of a set of data elements, each of the data elements modulating a carrier frequency of the signal, the data elements comprising pilots and so-called informative data elements. The invention is characterized in that said method comprises: a step which consists in determining a set of two-dimensional discrete prolate spheroidal sequences (DPSS); a step which consists in writing the transfer function in the form of a combination of at least some of the two-dimensional discrete prolate spheroidal sequences (DPSS); a step which consists in time-based and frequency-based two-dimensional interpolation of at least some of the coefficients of the combination, so as to obtain an estimation of the transfer function at any point of a time-frequency network.

[0001] The domain of this invention is transmission or broadcasting of digital data and/or sampled analogue data, particularly from and/or to mobile terminals.

[0002] Multicarrier modulations associated with error correction coding and interlacing have demonstrated their advantages particularly for high speed transmissions in a radio mobile environment as illustrated in patents FR-86 09622 and FR-95 05455 dealing with the OFDM (Orthogonal Frequency Division Multiplexing) modulation and IOTA (Isotropic Orthogonal Transform Algorithm) modulation.

[0003] In the case of a coherent demodulation of received data, it is essential to be able to make a good estimate of the channel at all points in a time-frequency plane, in order to correctly demodulate the information received by the receiver.

[0004] In particular, the invention relates to optimisation of channel estimating techniques in the case of a multiple carrier transmission, called multicarrier transmission.

[0005] Usually, when data are transmitted on a radio channel at high speed, the signal is subjected to Doppler effects (related to movements of the transmitter, receiver or any reflectors) and multiple path problems (for example related to reflection of the signal on different objects) causing a delay spread or variations of the signal amplitude. These effects depend particularly on the time and frequency considered. Thus, a channel with a variable time and frequency is complex and it has to be estimated correctly in order to be able to decode received data reliably.

[0006] According to the state of the art, in the case of a multicarrier modulation, the channel estimate from which degradations caused by a radio mobile channel can best be corrected, consists of the following steps:

[0007] insert reference carriers into the flow of information carriers at locations known to the receiver;

[0008] making use of values taken by these references in reception, to deduce the channel transfer function at these points in the time-frequency plane; and

[0009] starting from these results, obtain the channel transfer function on all points in a time-frequency network.

[0010] There are two particular methods of doing this:

[0011] the estimate by scattered pilots; in transmission, pilot symbols are distributed regularly in the time-frequency plane. An undersampled version of the channel is obtained starting from these pilots. In reception, a two-dimensional interpolation is carried out in time and frequency to determine the value of the channel at all points in the time-frequency network. This method is used particularly by the DVB-T (Digital Video broadcasting —Terrestrial) standard.

[0012] the estimate by reference multicarrier symbol (also called preamble); at least one reference symbol is placed at the beginning of a frame emitted on the channel. With this or these symbols, the channel is estimated in reception of the frame transmitted on each carrier. The channel can be considered as being quasi-static on a given frame (choice of system parameters such that the channel varies slowly compared with the symbol time), the channel estimate on the reference symbol(s) is valid for all OFDM symbols in the frame (in particular this method is applied to the ETSI (European Telecommunication Standard Institute) HIPERLAN/2 standard.

[0013] In the case of a scattered pilots estimate, the complexity of the channel estimator is usually limited by using two single-dimensional interpolation filters instead of a two-dimensional filter. The implementation complexity is then significantly lower for an acceptable degradation of the estimating quality.

[0014] For example for DVB-T, a Wiener filter is used for the frequency interpolation while the time interpolation is a simple linear interpolation. The linear interpolation is not optimum but its complexity is lower than the complexity of a two-dimensional FIR filter or two single-dimensional FIR filters in cascade.

[0015] Therefore, the interpolation is done in time and then in frequency.

[0016] Furthermore, it is difficult to make a Wiener filter adaptive, in other words to recalculate the coefficients of the reception filter as a function of the channel parameters. Consequently, it systematically becomes necessary to size the filter for the worst Doppler and delay spread case.

[0017] In particular, the purpose of the invention is to overcome these disadvantages of prior art.

[0018] More precisely, one purpose of the invention is to optimise decoding of a data frame emitted on a channel with a strong Doppler spread when multicarrier modulation is used.

[0019] Another objective is to optimise the estimate of the transfer function of a transmission channel within the framework of a multicarrier transmission.

[0020] Another purpose of the invention is to enable a fairly simple interpolation of the channel transfer function over an entire time-frequency network starting from estimated values of this function at locations of pilot symbols inserted in the transmitted data frame.

[0021] Another purpose of the invention is to offer a good compromise between the useful flow and the quality of the channel estimate.

[0022] These objectives, and others that will become clear later, are achieved by means of a method for estimating the transfer function for a multicarrier signal transmission channel formed by a series of symbols in time composed of a set of data elements, each data element modulating a carrier frequency of the signal, the data elements comprising firstly reference elements called pilots for which the transmission value is known to at least one receiver that will receive the signal, and secondly so-called information data elements representing at least one source signal to be transmitted.

[0023] According to the invention, this type of process comprises:

[0024] a step to determine a set of two-dimensional Discrete Prolate Spheroidal Sequences (DPSS);

[0025] a step to write the transfer function in the form of a combination of at least some of the two-dimensional discrete prolate spheroidal sequences in the set;

[0026] a two-dimensional interpolation step in time and in frequency of at least some coefficients of the combination, so as to obtain an estimate of the transfer function at all points in a time-frequency network.

[0027] Note that a time-frequency space is a limited space of model type used to represent a multi-carrier signal.

[0028] Thus, the invention enables the implementation of an estimate of a transfer function for a multicarrier signal transmission channel with reasonable complexity and that can be used to decode information data elements reliably and efficiently even in the presence of a large Doppler effect and/or delay spread of the channel.

[0029] According to one particular characteristic, the process is remarkable in that the write step uses a projection of pilots on two-dimensional discrete prolate spheroidal sequences.

[0030] Thus, the invention advantageously allows the use of a base composed of two-dimensional DPSSs.

[0031] According to one particular characteristic, during the determination step, the method uses at least one step for selection of at least some discrete prolate spheroidal sequences, as a function of at least one characteristic of the transmission channel.

[0032] Thus, the invention preferably enables the selection of two-dimensional DPSSs, particularly well adapted for the channel, which significantly improves its reliability.

[0033] In this manner, sizing of the filter for the worst Doppler and channel delay spread case is no longer systematically necessary since the estimate is adapted to the channel.

[0034] According to one particular characteristic, the process is remarkable in that the characteristic(s) of the transmission channel include the maximum Doppler frequency and/or the maximum channel delay spread.

[0035] Thus, the invention can advantageously take account of important characteristics of the channel present particularly for communications from and/or to mobile terminals.

[0036] According to one particular characteristic, the process is remarkable in that the selection step uses a sub-step to sort two-dimensional discrete prolate spheroidal sequences as a function of a predetermined energy criterion.

[0037] Thus, the invention advantageously only considers useful and/or necessary DPSSs. This limits the complexity of implementation while maintaining good reliability of the estimate of the transfer function.

[0038] According to one particular characteristic, the process is remarkable in that the number of two-dimensional discrete prolate spheroidal sequences selected during the selection step takes account of at least one transfer function estimate quality criterion.

[0039] Advantageously, the invention provides a means of estimating the transfer function in a manner suited to the required quality.

[0040] According to one particular characteristic, the process is remarkable in that the number of two-dimensional discrete prolate spheroidal sequences selected during the selection step is less than or equal to the number of signal pilots.

[0041] According to one particular characteristic, the process is remarkable in that during the determination step, a two-dimensional DPSS is obtained by taking a tensor product of at least two single-dimensional DPSSs.

[0042] Thus, the invention advantageously enables the use of two-dimensional DPSSs while limiting the complexity, by carrying out a judicious operation on at least two single-dimensional DPSSs. This is possible particularly if it is considered that the channel can be separated in time and in frequency.

[0043] According to one particular characteristic, the process is remarkable in that it uses an estimate-maximisation algorithm.

[0044] Thus, if the channel estimate quality obtained is not sufficient, the estimate-maximisation algorithm can advantageously improve it.

[0045] According to one particular characteristic, the process is remarkable in that during the writing step, it includes at least one step for estimating at least some coefficients of the combination using a least squares method.

[0046] The invention also uses a method for reception of a digital signal, remarkable in that it uses a step for estimating a signal transmission channel transfer function according to the method of estimating the transfer function of a multicarrier signal transmission channel.

[0047] The invention also relates to a multicarrier signal receiver formed from a time sequence of symbols composed of a set of data elements, each data element modulating a signal carrier frequency, the data elements comprising firstly reference elements called pilots for which the value on emission is known to at least one receiver that will receive the signal, and also information data elements representing at least one source signal to be transmitted.

[0048] The receiver is remarkable in that it includes:

[0049] means of memorising a set of two-dimensional discrete prolate spheroidal sequences (DPSSs);

[0050] means of writing a transfer function for a signal transmission channel in the form of a combination of at least some of the two-dimensional discrete prolate spheroidal sequences in the set;

[0051] means of making a two-dimensional interpolation in time and in frequency of at least some coefficients in the combination, so as to obtain an estimate of the transfer function at all points in the time-frequency network.

[0052] The invention also relates to an application of the method of estimating the transfer function of a multicarrier signal transmission channel in at least one of the following domains:

[0053] terrestrial digital broadcasting;

[0054] digital radio communications;

[0055] submarine data transmission.

[0056] The advantages of the reception process, the receiver and the application of the estimating process are the same as the advantages of the process for estimating the channel transfer function, and are not described in more detail.

[0057] Other characteristics and advantages of the invention will become more obvious after reading the following description of a preferred embodiment given as a simple illustrative and non-limitative example, and the attached drawings among which:

[0058]FIG. 1 illustrates two-dimensional DPSSs according to the invention;

[0059]FIG. 2 illustrates single-dimensional DPSSs used to make two-dimensional DPSSs as illustrated in FIG. 1;

[0060]FIG. 3 illustrates the fast decay of the eigenvalues of a single-dimensional DPSS in FIG. 2;

[0061]FIG. 4 shows a simplified block diagram of a transmitter station and a receiver station using two-dimensional DPSSs as illustrated with reference to FIG. 1;

[0062]FIG. 5 contains a flowchart used in a channel estimating module illustrated in FIG. 4;

[0063]FIG. 6 illustrates a radio communication model known per se;

[0064]FIG. 7 illustrates the result of the channel estimate used in FIG. 5 and applied to the model illustrated in FIG. 6.

[0065] The general principle of the invention is based on a two-dimensional interpolation technique concerning the scattered pilots estimate and more precisely the use of restrictions of two-dimensional Discrete Prolate Spheroidal Sequences (DPSS), the number of sequences being adapted to the characteristics of the channel considered and preferably not permanently fixed as a function of the worst propagation case.

[0066] The invention consists of using two-dimensional DPSSs that can be used in particular in cases in which the Doppler spread and Inter-Symbol Interference (ISI) of the channel (Doppler delay channel), such as satellite communication systems or mobile communication systems, and thus enable an optimum channel estimate in the case of a multicarrier system.

[0067] For a channel with a delay spread and Doppler effect, the channel transfer function that is locally stationary is modelled by narrow band white noise in two-dimensions (time and frequency). Its two-dimensional Doppler Power Spectrum is assumed to be flat (there is no prior information about the propagation channel), which encompasses all possible Doppler spread and time delay cases.

[0068] The size of the estimating window is limited and the number of two-dimensional DPSSs used is adapted as a function of the size of the estimating window.

[0069] Under these conditions, the use of two-dimensional DPSSs provides a means of obtaining a high quality channel estimate, while optimising the number of pilots for a reasonable implementation complexity. This estimating technique consists of interpolating the estimate of the coefficients of the discrete channel decomposition, obtained by pilot symbols based on two-dimensional DPSSs.

[0070] More precisely, in the two-dimensional case which is the purpose of this invention, the variable transfer function of the channel is modelled by two-dimensional narrow band noise. Since the receiver usually has incomplete knowledge of the characteristics of its local oscillator and of the delay-Doppler spectrum of the two-dimensional channel, this spectrum is assumed to be flat, with bounded support and with a width equal to the Doppler B_(D) spread on the frequency axis and the maximum delay T_(max) on the time axis. Therefore, it is modelled by a parallelepiped, defined as being the smallest parallelepiped encompassing the real support. Therefore, this model includes all possible delay and Doppler cases.

[0071] Self-correlation of T(f,t) for the Doppler power spectrum defined above, can be separated in time and frequency in the form: ${R_{T}\left( {{\Delta \quad f},{\Delta \quad t}} \right)}\overset{\Delta}{=}{E\left\lbrack {{{T\left( {f,t} \right)}T*\left( {{f + {\Delta \quad f}},{t + {\Delta \quad t}}} \right\rbrack}\quad = {N_{c}{{sinc}\left( {\pi \quad \Delta \quad f\quad T_{\max}} \right)}{{sinc}\left( {2\pi \quad \Delta \quad t\quad f_{D}} \right)}}} \right.}$

[0072] where T_(max) is the maximum spread, 2f_(D) is the Doppler spread, f is a frequency, Δf is a frequency variation, t is a time, Δt is a time variation and sinc is the cardinal sine function.

[0073] We will consider a block with dimensions N₁×N₂, where N₁ (number of time symbols, for example equal to 25) corresponds to the time dimension and N₂ (number of sub-carriers, for example equal to 35) corresponds to the frequency dimension. The parameters N₁ and N₂ are predetermined as a function of a given complexity and/or a compromise between the complexity and the required reliability of the estimate. The DPSSs are two-dimensional functions defined on the product space C^(N1)×C^(N2), where C^(N1) corresponds to the complex time space with N₁ dimensions and C^(N2) corresponds to the frequency space with N₂ dimensions.

[0074] Due to the separable nature of the channel in time and in frequency, these two-dimensional DPSSs are defined as the tensor product of single-dimensional DPSSs.

[0075] Thus defined, 2D (in other words two-dimensional) DPSSs can be used to obtain decorrelated channel decomposition coefficients:

[0076] Let R_(c) be the covariance vector of the channel, defined by: $\underset{\underset{\_}{\_}}{\underset{\underset{\_}{\_}}{R_{c}}}\quad \underset{\underset{\_}{\_}}{\Delta}\left\{ {E\left\lbrack {c_{ij}c_{kl}^{*}} \right\rbrack} \right\}_{jkl}$

[0077] where:

[0078] i and k are taken from 1 to N₁ and j and 1 are taken from 1 to N₂;

[0079] C_(k1)* represents the conjugate of the complex value C_(k1);

[0080] the function E[x] representing the expectancy or the average of the variable x; and

[0081] the complex vector c to be estimated in this case representing an implementation of the two-dimensional discrete channel at the emitted data stream, and being composed of N₁.N₂ complex coefficients c_(ij).

[0082] The channel can be separated, which is equivalent to writing that the vector R_(c) is the tensor product of covariance matrices corresponding to each of the time and frequency dimensions: $\begin{matrix} {\underset{\underset{\_}{\_}}{\underset{\underset{\_}{\_}}{R_{c}}} = \left\{ {N_{c}\quad \sin \quad {c\left( {{\pi \left( {j - 1} \right)}v_{0}T_{\max}} \right)}\sin \quad {c\left( {2\quad {\pi \left( {i - k} \right)}{\tau_{0} \cdot f_{D}}} \right)}} \right\}_{ijkl}} \\ {= {N_{c}\quad {\underset{\underset{\_}{\_}}{R^{(1)}} \otimes \quad \underset{\underset{\_}{\_}}{R^{(2)}}}}} \end{matrix}$ where $\underset{\underset{\_}{\_}}{R^{(1)}} = \left\{ {{\sin \quad {c\left( {{\pi \left( {j - l} \right)}{v_{0} \cdot T_{\max}}} \right\}}_{jl}};{{{and}\quad {where}\underset{\underset{\_}{\_}}{R^{(2)}}} = \left\{ {\sin \quad {c\left( {2\quad {\pi \left( {i - k} \right)}{\tau_{0} \cdot f_{D}}} \right\}}_{ik}} \right.}} \right.$

[0083] Therefore, the decomposition coefficients of the discrete channel based on 2D DPSSs are decorrelated.

[0084] It can be demonstrated that the 2D DPSSs thus created are actually eigenvectors of the order 4 covariance vector of the two-dimensional channel.

[0085] It is possible to find an orthonormed base for the complex space with N₁.N₂ dimensions C^(N1.N2), for which the normalized vectors {p ^((i,j))}, i=l . . . N₁ and j=l . . . N₂ (where the notation x=a . . . b means that values of x are taken from a to b) are the eigenvectors of the covariance vector.

[0086] Therefore the family of two-dimensional DPSSs is the family $\left\{ \left( {{\underset{\underset{\_}{\_}}{p}}^{({i,j})},{\lambda^{({l,i})}\lambda^{({2,j})}}} \right) \right\}_{i = {{1\quad j} = 1}}^{N_{1}N_{2}}$

[0087] where:

[0088] the orthonormed two-dimensional DPSSs are the following matrices ${\underset{\underset{\_}{\_}}{p}}^{({i,j})} = \left\{ \left( {{Pk}^{({1,i})},{Pl}^{({2,j})}} \right) \right\}_{k = {{1\quad l} = 1}}^{N_{1}N_{2}}$

[0089]  and the associated eigenvalues are λ^((1,i)).λ^((2,j)) where: $\left\{ \left( {{\underset{\_}{p}}^{({1,i})},\lambda^{({1,i})}} \right) \right\}_{i = 1}^{N_{1}}\quad {and}\quad \left\{ \left( {{\underset{\_}{p}}^{({2,i})},\lambda^{({1,i})}} \right) \right\}_{j = 1}^{N_{2}}$

[0090] two families of single-dimensional DPSSs on the time dimension and the frequency dimension respectively, and their associated eigenvalues.

[0091] The $\left\{ {\underset{\_}{p}}^{({1,i})} \right\}_{i = 1}^{N_{1}}$

[0092] sequences are time sequences with size N₁ with a limited normalized band [−W_(p) W_(t)] and the sequences $\left\{ {\underset{\_}{p}}^{({2,j})} \right\}_{j = 1}^{N_{2}}$

[0093] are frequency sequences with size N₂ with normalized limited band [0; 2W_(f)] where W_(t)=f_(D)×τ₀ and W_(f)= $\frac{T_{\max} \times v_{0}}{2}.$

[0094] Vectors corresponding to the time dimension are equal to the restriction of N₁ DPSSs, real discrete time sequences with dimension N₁ and with a normalized limited band [−W_(p); W_(t)] concentrated mostly on the discrete interval ${\overset{i}{n}}_{n = 0}^{N_{1} - 1}$

[0095] corresponding to the restriction to the time support of the transmitted block. By definition, the most concentrated DPSSs are those with the greatest energies within the interval considered.

[0096] Similarly, vectors corresponding to the frequency dimension are equal to the restriction of the N₂ DPSSs, real discrete frequency sequences with dimension N₂ and with a limited normalized band [0; 2W_(f)] with their greatest concentration on the discrete interval ${\overset{i}{n}}_{n = 0}^{N_{2} - 1}$

[0097] corresponding to the restriction to the frequency support of the transmitted block.

[0098] The decomposition coefficients of the two-dimensional discrete channel are then given by: $\underset{\underset{\_}{\_}}{c} = {\sum\limits_{i = 1}^{N_{1}}{\sum\limits_{j = 1}^{N_{2}}{b_{ij}{\underset{\underset{\_}{\_}}{p}}^{({i,j})}}}}$

[0099] For reasons of clarity, in the remaining part of this description we will use a single index giving: $\underset{\underset{\_}{\_}}{c} = {\sum\limits_{i}{b_{i}{\underset{\underset{\_}{\_}}{p}}^{(i)}}}$

[0100] Therefore, the choice of the right single-dimensional functions provides a means of building good two-dimensional interpolation functions.

[0101] Only the most concentrated DPSSs are necessary to describe the channel. Only the most concentrated N′ 2D DPSSs are kept to interpolate the channel.

[0102] According to the preferred embodiment, the N′ most concentrated 2D DPSSs are obtained by calculating the N₁ eigenvalues λ^((1,j)) in the time dimension, and the N₂ eigenvalues λ^((2,i)) in the frequency dimension, and then selecting the N′ highest values of the products λij=λ^((1,j)).λ^((2,j)).

[0103] According to a variant that is less complex but gives a worst estimate (for a given number of selected 2D DPSSs) than the preferred embodiment, the selection is simplified by firstly selecting the N′₁ highest values for λ^((1,j)) and the N′₂ highest values for λ^((2,j)), to obtain N′ (in this case equal to N′₁.N′₂) values giving the most concentrated 2D DPSSs.

[0104] Therefore, the channel estimate corresponds to:

ĉ ^((N′))=c+n _(c)

[0105] where n _(c) is the estimate noise.

[0106] By analogy with the signal to noise ratio, the quality criterion of the following estimate is defined as follows: ${Q\left( N^{\prime} \right)} = {\frac{E\left\lbrack {\underset{\underset{\_}{\_}}{c}}^{2} \right\rbrack}{E\left\lbrack {{\underset{\underset{\_}{\_}}{n}}_{c}}^{2} \right\rbrack}\quad {where}\quad {.}}$

[0107] defines the norm of a matrix.

[0108] Therefore we can write: ${Q\left( N^{\prime} \right)} = \frac{\sum{{Eigenvalues}\quad {of}\quad {all}\quad {SSADs}}}{\sum{{Eigenvalues}\quad {of}\quad {rejected}\quad {SSADs}}}$

[0109] In the first step of the estimating algorithm, the coefficients are obtained by projecting reference symbols on DPSSs that are kept. Since not all DPSSs are kept to estimate the channel, these projections give an acceptable approximation of the channel decomposition.

[0110] In one preferred embodiment, we restrict the discussion to the projection of reference symbols on kept DPSSs without using an EM (Estimate—Maximisation) algorithm described later.

[0111] Nevertheless, as a variant, in order to optimise the channel estimate for arbitrary synchronization reasons and taking account of all or some of the coded structure of the transmitted data, an estimate—maximisation algorithm is used in a manner similar to that presented in patent FR-96 05200 (deposited by the same applicants and entitled “Signal Numérique à blocs de reference multiples pour l'estimation de canal, procédés d'estimation de canal et récepteurs correspondants” (“Digital signal with multiple reference blocks for channel estimating, channel estimating methods and corresponding receivers”). In this way, not only are the reference sub-carriers used, but useful carriers are also used in order to improve the estimator quality in terms of the Maximum A Posteriori (MAP) probability criterion.

[0112] As an illustration, FIG. 1 shows DPSSs 15 to 18 on graphs 11 to 14. For example, in this case we will consider the case of a “Vehicular B at 250 km/h” type model where f_(D) is close to 930 Hz, T_(max) equals 20 μs, τ₀ is equal to 133.33 μs and v₀ is equal to 3.75 kHz. It is deduced that time is limited to one interval [−Wt, +Wt] where Wt is equal to 0.0622 and frequencies are limited to an interval [0, +Wf] where Wf is equal to 0.0375. According to one example considered, in this case it was decided to represent the four two-dimensional DPSSs with the highest energy (the energy of one DPSS being equal to the product of the eigenvalues of the two corresponding 1D DPSSs, namely λ^((1,i)) λ(^(2,j))) depending on their amplitude (vertical axis) as a function of the time and frequency (horizontal axes) (whereas in fact, for example, the 25 highest energy two-dimensional DPSSs were kept to estimate the channel).

[0113] We will now describe details of how single-dimensional DPSSs families are obtained.

[0114] Consider the case of the time dimension (the frequency dimension is similar).

[0115] The complex vector c corresponding to the time dimension to be estimated in this case represents an embodiment of the discrete channel at the transmitted data stream, and in this case is composed of N (equal to N₁) complex coefficients c_(n).

[0116] It is possible to find an orthonormed base for the complex space with N dimensions C^(N), for which the normalized vectors {p^((i))}, i=1 . . . N are the eigenvectors of the covariance matrix of c associated with the eigenvalues {λ_(i)}, i=1 . . . N.

[0117] The Doppler spectrum of the channel in this case is assumed to be flat with a bounded support and a width equal to the Doppler spread B_(D)=2f_(D), f_(D) representing the maximum Doppler frequency.

[0118] In this case, these vectors are equal to the restriction of N DPSSs, real discrete sequences with dimension N and with a normalized limited band [−W′;W′] where W′=W_(t)=the most concentrated f_(D)×ρ₀ in the discrete interval {n}_(n = 0)^(N − 1),

[0119] corresponding to the restriction to the time support of the transmitted block. By definition, the most concentrated DPSSs are the DPSSs with the highest energies within the interval considered. Thus, we only keep the N DPSSs with the highest energy and the others are rejected.

[0120] In other words, $\underset{\_}{c} = {\sum\limits_{i = 1}^{N}{b_{i}{\underset{\_}{p}}^{(i)}}}$

[0121] where the factors bi are independent complex gaussian random variables satisfying the following relations:

[0122] the variance of the variables b_(i) is equal to λ₂:

E[|b_(i)|²]=λ_(i)  (relation(1))

[0123] the covariance between b_(i) and b_(j), where i and j are different is zero:

E[b_(i)b_(j)*]=0 if i≠j

[0124] where the variable b_(j)* is the conjugate of the complex variable b_(j).

[0125] The basic vectors (in other words the DPSSs that were kept) are given by the relation:

E[c.c ^(H)].p^((i))=λ_(i)p^((i))

[0126] where (.)^(H) is the transposition—conjugation operator of the vector (.).

[0127] The {p^((i))}, i=1 . . . N vectors satisfy the following system of equations: ${2W^{\prime}{\sum\limits_{j = 1}^{N}{\sin \quad {c\left( {2{W^{\prime}\left( {n - j} \right)}} \right)}p_{j}^{(i)}}}} = {\lambda_{i}{\underset{\_}{p}}_{n}^{(i)}}$ n = 1  …  N

[0128] which corresponds to the definition of single-dimensional DPSSs with size N and a normalized band [−W′;W′].

[0129] For example, FIG. 2 shows an example of single-dimensional DPSSs. The 25 highest energy DPSSs are selected in this example (N is equal to 25). Graph 20 shows the three DPSSs 21, 22 and 23 with the highest energy for which the amplitude (ordinate axis) is given as a function of time (abscissa axis for which the scale is in normalized units, in other words as a number of symbol times). The concentration or energy of the DPSSs depends on the eigenvalues; as the modulus of an eigenvalue corresponding to a DPSS increases, the energy of the corresponding DPSS also increases. Note that the three DPSSs 21, 22 and 23 sorted in decreasing order of energy correspond to three eigenvalues equal to 0.99924, 0.97646 and 0.77680 respectively (which are therefore also sorted in decreasing order).

[0130] Since the concentration of a DPSS is equal to the n^(th) eigenvalue λ_(n) of the covariance matrix of c, relation (1) shows that only the most concentrated DPSSs are necessary to describe the channel.

[0131] It is assumed that these eigenvalues are stored in decreasing order. The article “Prolate Spheroidal Wave Function, Fourier Analysis and Uncertainty-I” written by D. SLEPIAN and H. O. POLLAK published in the “The Bell System Technical Journal” in January 1961, demonstrates that these eigenvalues have a very fast decay and quickly disappear after 2W′N values. Note that according to the example illustrated in FIG. 3, the product 2W′N is equal to 3.1 for values of W′ and N equal to 0.062 respectively, and that the eigenvalues expressed in dB (in other words according to a 10 Log(λ²/λ₀ ²) type scale, where λ represents the value considered and λ₀ is the largest eigenvalue) decay very quickly after this threshold of 3.1.

[0132] The DPSS family is orthogonal. However, there is nothing to assure that undersampled DPSSs are orthogonal. Therefore, projection of reference symbols onto these DPSSs gives an approximation of the channel decomposition.

[0133] The following three steps are carried out during the channel estimate:

[0134] the receiver uses N′ reference symbols of the modulated block to calculate the coefficient di corresponding to the channel decomposition on the reduced family of N′ DPSSs. Due to the presence of an additive noise, the following weighting coefficients w_(i) are applied to the result obtained: $w_{i} = \frac{\lambda_{i}}{\lambda_{i} + \frac{N_{0}}{N_{c}E_{s}}}$

[0135] where N_(c) is the channel variance, N₀ is the additive noise variance, and E_(s) is the energy of the transmitted reference symbols.

[0136] The result for the estimate of the decomposition coefficients is then:

[0137] {circumflex over (b)}_(i)=w_(i)d_(i) i=1 . . . N′

[0138] The EM (Estimate —Maximisation) algorithm is used to improve the estimates obtained. Independence of the coefficients bi simplifies calculations of the maximisation step. To improve the estimates obtained, the EM algorithm described in detail in patent FR-96 05200 mentioned above uses not only the reference sub-carriers but also the useful carriers in order to improve the estimator quality in terms of the Maximum A Posteriori (MAP) probability criterion.

[0139] The receiver then determines a channel estimate by making an interpolation using the kept DPSSs.

[0140] The reader can refer to patent FR-96 05200 mentioned above for further details about use of single-dimensional DPSSs.

[0141] We will now present an embodiment of the channel estimate according to the invention, with reference to FIG. 4.

[0142]FIG. 4 shows a simplified diagram illustrating particularly:

[0143] a transmitting station 41 transmitting a signal formed of a modulated data stream according to the invention;

[0144] a transmission channel 42 for the transmitted signal; and

[0145] a receiving station 43 for the transmitted signal, particularly including a channel estimate according to the invention.

[0146] In particular, the transmitting station 41 comprises an information source 411 with an arbitrary flow generating binary or non-binary data corresponding to an arbitrary type of source signals (sound, picture, data, etc.). These data may be subjected to source and channel coding 412. Channel coding uses an error correction code adapted to the transmission channel 42.

[0147] Coded data generated from these codes (useful symbols) are then organized into data streams, and are modulated (413). Therefore, they are suitably distributed and interleaved onto several data streams in order to obtain the necessary diversity and to decorrelate fading affecting the transmitted symbols. Reference elements are also introduced into each data stream, according to distribution principles specified later. Finally, the data are modulated according to OFDM/QAM (Quadrature Amplitude Modulation) type modulation (for example OFDM with guard interval), or OFDM/OQAM (Offset Quadrature Amplitude Modulation) (for example the IOTA modulation). Note that reference symbols are inserted in the block to be transmitted after the “mapping” and before the OFDM modulator. These symbols are uniformly distributed in the time-frequency plane. Refer to patents FR-86 09622 and FR-95 05455 for a detailed description of the modulation system.

[0148] The signal is then transposed in frequency, and is amplified and transmitted by the transmitter 414 through the channel 42.

[0149] The channel 42 is an arbitrary radio transmission channel. The signal passes in channel 42 and in this channel, for example, noise may be added, it may follow multiple paths, or interference may be introduced, or it may be influenced by the Doppler effect.

[0150] Those skilled in the wireless transmission field usually models the channel, particularly so as to optimise coding/decoding and modulation/demodulation.

[0151] In this case, according to one model of this channel, it is assumed that the channel is affected by:

[0152] white gaussian noise characterised by a Signal to Noise Ratio (SNR) assumed to be known;

[0153] the Doppler effect, for which the power spectrum is normalized by a parallelepiped of normalized half bands Wt and Wf and

[0154] delay spread (particularly related to multiple paths resulting from signal reflections on environmental elements (for example buildings)).

[0155] Some of these models are normalized. For example, it can be assumed that channel 42 can be modelled using the “Vehicular B” type at 250 km/h described in the ETSI report “TR 101 112 V3.2.0 appendix B.1 ETSI” that is typical of a channel between a fixed station and a vehicle moving at a speed equal to 250 km/h.

[0156] Input stages in the receiving station 43 are conventional. The receiving station 43 receives the signal sent by the transmitting station 41 and transmits it through the channel 42. The signal corresponding to a received data stream is preamplified 431 and is then converted into an intermediate frequency in order to achieve tuned filtering of channel 432. The intermediate frequency signal is then converted into base band on two channels in quadrature, and is then sampled (437).

[0157] The receiver uses the samples corresponding to a data stream to determine an estimate 438 described in detail below.

[0158] This estimate enables a reliable demodulation 439 of the samples.

[0159] The demodulated data are then decoded, possibly de-interleaved to be transmitted to the recipient 434.

[0160] The channel estimate also enables automatic gain correction (AGC) 435 that controls the preamplification 431.

[0161] In this case, it is assumed that the data streams are short. There is no need for synchronization. As a variant, for long data streams, synchronization taking account of the samples makes it possible to synchronise correctly on transmitted symbols and/or the data stream itself to sample the received signal correctly.

[0162] We will now describe the channel estimate 438 according to the invention in detail, with reference to the flowchart described in FIG. 5.

[0163] The two-dimensional DPSSs adapted to the channel type used are determined in an initialisation step 51, using the following method.

[0164] It is assumed that system parameters (in other words the symbol time τ₀ and the space between subcarriers ν₀) were chosen such that the channel 42 can be considered as being quasi-constant on the scale of an OFDM mesh (τ₀, ν₀) (an OFDM mesh corresponds to an area with dimension τ₀ in time and ν₀ in frequency in the time-frequency network), namely: $\quad\left\{ \begin{matrix} {{v_{0}{\operatorname{<<}\frac{1}{T_{Max}}}} = B_{c}} \\ {{\tau_{0}{\operatorname{<<}\frac{1}{2f_{D}}}} = T_{c}} \end{matrix} \right.$

[0165] where

[0166] Bc represents the channel consistency band; and

[0167] Tc represents the channel consistency time.

[0168] (the “<<” symbol means “very much less than”).

[0169] The channel then behaves like a multiplication channel, characterised by an amplitude and a phase corresponding to the value of T(f,t), a variable transfer function of the channel, for the time and frequency considered.

[0170] Therefore, for the mesh corresponding to a transmitted symbol α_(m,n) we have:

[0171] T(f,t)=T(mν₀, nτ₀) and

[0172] τ(mν₀,nτ₀)=ρ_(m,n)e^(iθ) ^(_(m,n))

[0173]  where ρ represents an amplitude and θ represents a phase.

[0174] Let us use the following matrix notation:

[0175]α ^(Δ) ={α_(m,n)}_(m,n) transmitted binary symbols matrix

[0176]r ^(Δ) ={r_(m,n)}_(m,n) matrix of samples received at the output from the demodulator

[0177]c ^(Δ) ={T(mν₀, n_(τ) ₀)}_(m,n) matrix of discrete channel coefficients

[0178]n ^(Δ) ={n_(m,n)}_(m,n) matrix of coefficients of discrete complex gaussian additive white noise

[0179] The discrete channel is modelled by the following relation:

r=c·α+n

[0180] where the “·” operator represents the term-by-term product.

[0181] Similarly, since the channel is locally considered as a multiplication channel:

r _(m,n) =C _(m,n)α_(m,n) +n _(m,n)

[0182] It is required to determine the values C_(m,n).

[0183] Remember that the reference symbols are inserted in the block to be transmitted before the OFDM modulator and these symbols are uniformly distributed in the time-frequency plane.

[0184] The characteristics of the propagation channel, namely maximum Doppler spread (bilateral) and maximum delay spread (unilateral) are assumed to be known (for example, these values can be estimated using a channel model and/or they can be measured).

[0185] It is assumed that the Doppler Power Spectrum of the channel is modelled by a parallelepiped of normalized half-bands W_(t) and W_(f) (as specified in the description of the general principle of the invention):

W _(t) =f _(D) ×τ ₀

[0186] and $W_{f} = \frac{T_{\max} \times v_{0}}{2}$

[0187] The single-dimensional DPSSs corresponding to these normalized bands are calculated in advance:

[0188] In the time dimension:

[0189] let R (N, W_(t)) be the covariance matrix with dimension N (corresponding to the time dimension of the treated block) of complex white noise with 2W_(t) normalized bands,

R _(m,n=)2W_(t) sin (2πW _(t)(m−n))

[0190] 1D DPSSs (in other words single-dimensional) corresponding to the time dimension are defined as being normalized eigenvectors of the matrix R(N,W_(t)), sorted in increasing order of their eigenvalues.

[0191] In the frequency dimension:

[0192] The procedure is similar:

[0193] Let R (N, Wf) be the covariance matrix with dimension N (corresponding to the frequency dimension of the block being considered) of complex white noise with 2Wf normalized bands,

R_(m,n)=2W_(f) sin (2πW _(f)(m−n))

[0194] 1D DPSSs corresponding to the frequency dimension are defined as being normalized eigenvectors of the matrix R(N,W_(f)) sorted in increasing order of their eigenvalues.

[0195] The 2D (in other words two-dimensional) DPSSs are then constructed using term by term multiplication of 1D DPSSs corresponding to a time dimension and a 1D DPSS corresponding to a frequency dimension:

[0196] P_(ij)=P_(i)P_(j) ^(T); and

[0197] λ_(ij)=λ_(i)λ_(j where:)

[0198] P_(ij) defines a two-dimensional DPSS of eigenvalues

[0199] P_(i) and P_(j) correspond to single-dimensional DPSSs with the corresponding eigenvalues λ_(i) and λ_(j) (notations simplified to facilitate reading).

[0200] The DPSSs with the highest energies among the resulting DPSSs are selected.

[0201] According to a first embodiment, the number N′ of 2D DPSSs selected is determined in advance and may for example be equal to approximately 25.

[0202] According to a second embodiment, the number N′ of selected 2D DPSSs is such that the quality factor defined as being the ratio of the sum of eigenvalues to the sum of rejected eigenvalues is greater than a given quality factor Q₀ (for example of the order of 30 dB). In summary, the number N′ is determined by the following criterion: $Q = {\frac{\sum{{all}\quad {eigenvalues}}}{\sum{{rejected}\quad {eigenvalues}}} > Q_{0}}$

[0203] Thus, for a time-frequency block with a given size and given characteristics of the propagation channel, the number of two-dimensional DPSSs necessary for interpolation of the complex gain of the channel is generated from the estimate of this complex gain obtained at the locations of the pilot symbols.

[0204] Thus, during the initialisation step 51, the selected two-dimensional DPSS parameters are stored in the channel estimating module 438 of the receiving station.

[0205] Note that the two-dimensional DPSSs may be retained during the design of the channel estimating module 438 and may be used statically during a channel estimate.

[0206] According to a more complex variant (the complexity remaining reasonable, due to the method of building up 2D DPSSs obtained by the tensor product of two 1D DPSSs), the 2D DPSSs are determined dynamically by the channel estimating module 438.

[0207] According to this variant, the 2D DPSSs to be used are chosen as a function of the characteristics of the propagation channel (particularly the maximum propagation time measured starting from the time elapsed between reception of the first symbol and the last echo corresponding to this symbol).

[0208] Then, during a step 52, the receiving station 43 waits for and then receives data blocks to be demodulated. In particular, this block includes pilots that will be used for the channel estimate.

[0209] Then, pilots are projected onto the selected two-dimensional DPSSs during a combination write step 53.

[0210] In order to make this projection, the first step is to calculate the channel decomposition coefficients (d_(i))_(i=1, . . . , K) in the space containing 2D DPSSs kept starting from the transmitted reference symbols (a_(i))_(i=1, . . . , K): $\left\lceil \begin{matrix} d_{1} \\ \vdots \\ d_{N^{\prime}} \end{matrix} \right\rceil = {\left\lbrack {P_{1}\quad \cdots \quad P_{N^{\prime}}} \right\rbrack^{\#}\left\lceil \begin{matrix} a_{1} \\ \vdots \\ a_{k} \end{matrix} \right\rceil}$

[0211] where:

[0212] N′ is the number of selected DPSSs;

[0213] K is the number of pilots;

[0214] p is a 2D DPSS (N′ 2D DPSSs kept) and p# is its pseudo inverse, as follows: $\left\lceil \begin{matrix} a_{1} \\ \vdots \\ a_{k} \end{matrix} \right\rceil = {\left\lbrack {P_{1}\quad \cdots \quad P_{N^{\prime}}} \right\rbrack \cdot \left\lceil \begin{matrix} d_{1} \\ \vdots \\ d_{N^{\prime}} \end{matrix} \right\rceil}$

[0215] The number K of pilot symbols used must be greater than or equal to the number N′ of DPSSs kept. Remember that P_(i(1≦i≦N′)) is a 2D DPSS, namely a vector with dimension K. Therefore, since K is greater than or equal to N′, the number of channel decomposition coefficients d_(i) is the same as the number of 2D DPSSs.

[0216] Pseudo-inverses are stored in a non-volatile memory (ROM) on the receiving station, for each processed block size (an estimating window size is usually defined corresponding to the size of the smallest block that can be transmitted) and each propagation environment.

[0217] Since the channel is noisy and the SNR (Signal/Noise Ratio) is assumed to be known, the result is an estimate of the channel decomposition coefficients (a_(i))_(i=1, . . . , K) in the 2D DPSSs space using the least squares method (LSM). ${\hat{a}i} = {\frac{\lambda_{i}}{\lambda_{i} + \frac{N_{0}}{E_{s}}}d_{i}}$

[0218] λ_(i) is the eigenvalue of the ith 2D DPSS;

[0219] N₀ is the noise variance; and

[0220] E_(s) is the energy of pilots on reception.

[0221] The channel estimate Ĉ is then given by the following equation: $\hat{C} = {\left\lbrack {P_{1}\quad \cdots \quad P_{N^{\prime}}} \right\rbrack^{\#}\left\lceil \begin{matrix} {\hat{a}}_{1} \\ \quad \\ {\hat{a}}_{k} \end{matrix} \right\rceil}$

[0222] Then, during a step 54, the channel estimating module 438 makes an estimate of the channel transfer function by interpolation of the estimate made on the pilots. Pilot symbols known to the receiver are used to obtain a discrete channel decomposition based on the 2D DPSSs. Knowledge of the channel at any point in the time-frequency network can then be obtained by transposition in the time base.

[0223] As an illustration, FIG. 6 show two-dimensional signals for a “Vehicular B” type channel at 250 km/h defined by the ETSI under reference “TR 101 112 V.3.2.0 appendix B.1 ETSI”.

[0224] The precise channel characteristics are as follows:—

[0225] The channel is a rural channel

[0226] The Doppler spectrum is a horn (spectrum corresponding to a power spectrum of a pure sine curve after propagation in a mobile channel with reflections assumed to be isotropic).

[0227] Delays (due to multiple paths) are summarised in the following table: Delay in nsec Attenuation in dB 0 −2.5 300 0 8900 −12.8 12900 −10.0 17100 −25.2 20000 −16.0

[0228] The modulation is of the OFDM/IOTA type with the following parameters:

[0229] the IOTA symbol type is equal to 133.33 μs; and

[0230] the space between subcarriers is equal to 3.75 kHz.

[0231] It is assumed that the channel is not noisy.

[0232]FIG. 7 illustrates the same channel as that illustrated in FIG. 6 but estimated according to the invention using 2D DPSSs.

[0233] The characteristics of 2D DPSSs are as follows:

[0234] the number of 2D DPSSs kept is equal to 25

[0235] the time pitch of pilot symbols is equal to 5 (so that four useful symbols can be obtained for one time reference symbol)

[0236] the frequency pitch of pilot symbols is equal to 8 (so that seven useful symbols can be obtained for one frequency reference symbol).

[0237] It can be seen that the channel according to the model described in FIG. 6 and its estimate according to invention are very similar. Thus, the channel estimate with two-dimensional DPSSs is very reliable and efficient.

[0238] Obviously, the invention is not restricted to the example embodiments mentioned above.

[0239] In particular, those skilled in the art could add any variant to the type of multicarrier modulation, particularly OFDM, which for example could be of the OFDM/QUAM or OFDM/OQAM type.

[0240] Note also that the invention can be used for applications in a wide variety of domains, particularly when high spectral efficiency is required and the channel is highly non-stationary, particularly when it is found that it would be judicious to use multicarrier modulation techniques.

[0241] A first category of applications relates to terrestrial broadcasting (for example, Digital Audio Broadcasting (DAB)), particularly for pictures, sound and/or data.

[0242] A second application category relates to digital radiocommunications. In particular, the invention can be used for applications in high-speed digital communication systems originating from or leading to mobiles (for example according to third generation mobile communication standards) and high speed LANs using multicarrier modulation techniques.

[0243] A third category of applications is submarine transmissions for which multicarrier modulation techniques are very suitable.

[0244] Therefore, in general, the invention can be used for applications in all domains in which multicarrier modulation techniques themselves are used (for example in systems combining Code Division Multiple Access (CDMA) and QFDM (particularly (multicarrier CDMA (MC-CDMA)). 

1. Method for estimating the transfer function for a multicarrier signal transmission channel formed by a series of symbols in time composed of a set of data elements, each of the data elements modulating a carrier frequency of the signal, the data elements comprising reference elements for which the transmission value is known to at least one receiver that will receive the signal, and information data elements representing at least one source signal to be transmitted, characterized in that it comprises: a step to determine a set of two-dimensional Discrete Prolate Spheroidal Sequences (DPSS); a step to write the transfer function in the form of a combination of at least some of the two-dimensional discrete prolate spheroidal sequences in the set; a two-dimensional interpolation step (54) in time and in frequency of at least some coefficients of the said combination, so as to obtain an estimate of the said transfer function at all points in a time-frequency network.
 2. Method according to claim 1, characterized in that the said write step uses a projection of the reference elements on the two-dimensional discrete prolate spheroidal sequences.
 3. Method according to claim 1, characterized in that during the determination step, the method uses at least one step for selection of at least some of the discrete prolate spheroidal sequences (DPSSs), as a function of at least one characteristic of the transmission channel.
 4. Method according to claim 3, characterized in that the characteristic(s) of the transmission channel include the maximum Doppler frequency (f_(d)) and/or the maximum channel delay spread (T_(max)).
 5. Method according to claim 3, characterized in that the selection step uses a sub-step to sort the two-dimensional discrete prolate spheroidal sequences (DPSSs) as a function of a predetermined energy criterion (Q₀).
 6. Method according to claim 3, characterized in that the number of two-dimensional discrete prolate spheroidal sequences selected during the selection step takes account of at least one estimate quality criterion of the transfer function.
 7. Method according to claim 3, characterized in that the number of two-dimensional discrete prolate spheroidal sequences selected (N′) during the selection step is less than or equal to the number of reference elements (K) of the signal.
 8. Method according to claim 1, characterized in that during the determination step, a two-dimensional DPSS of the set is obtained by taking a tensor product of at least two single-dimensional DPSSs.
 9. Method according to claim 1, characterized in that it uses an estimate-maximisation (EM) algorithm.
 10. Method according to claim 1, characterized in that during the said writing step, it includes at least one step for estimating at least some coefficients of the combination using a least squares method.
 11. Method for reception of a digital signal, characterized in that it uses a step for estimating a transfer function of a transmission channel of the said signal using the method according to claim
 1. 12. Multicarrier signal receiver formed from a time sequence of symbols composed of a set of data elements, each of the data element modulating a carrier frequency of the signal, the data elements comprising reference elements for which the value on emission is known to at least one receiver that will receive the signal, and also information data elements representing at least one source signal to be transmitted, characterized in that it includes: means of memorizing a set of two-dimensional discrete prolate spheroidal sequences (DPSSs); means of writing a transfer function for a transmission channel of the said signal in the form of a combination of at least some of the two-dimensional discrete prolate spheroidal sequences in the set; means of making a two-dimensional interpolation in time and in frequency of at least some coefficients in the combination, so as to obtain an estimate of the said transfer function at all points in the time-frequency network.
 13. Application of the estimating method according to claim 1 to at least one of the following domains: terrestrial digital broadcasting; digital radiocommunications; submarine data transmission.
 14. Method according to claim 2, characterized in that during the determination step, the method uses at least one step for selection of at least some of the discrete prolate spheroidal sequences (DPSSs), as a function of at least one characteristic of the transmission channel.
 15. Method according to claim 14, characterized in that the characteristic(s) of the transmission channel include the maximum Doppler frequency (f_(d)) and/or the maximum channel delay spread (T_(max)).
 16. Method according to claim 15, characterized in that the selection step uses a sub-step to sort the two-dimensional discrete prolate spheroidal sequences (DPSSs) as a function of a predetermined energy criterion (Q₀).
 17. Method according to claim 14, characterized in that the selection step uses a sub-step to sort the two-dimensional discrete prolate spheroidal sequences (DPSSs) as a function of a predetermined energy criterion (Q₀).
 18. Method according to claim 14, characterized in that the number of two-dimensional discrete prolate spheroidal sequences selected during the selection step takes account of at least one estimate quality criterion of the transfer function.
 19. Method according to claim 15, characterized in that the number of two-dimensional discrete prolate spheroidal sequences selected during the selection step takes account of at least one estimate quality criterion of the transfer function.
 20. Method according to claim 14, characterized in that the number of two-dimensional discrete prolate spheroidal sequences selected (N′) during the selection step is less than or equal to the number of reference elements (K) of the signal. 